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ABSTRACT 

We use mid-infrared spectroscopy of unobscured active galactic nuclei (AGNs) 
to reveal their native dusty environments. We concentrate on Seyfert 1 galax- 
ies, observing a sample of 31 with the Infrared Spectrograph aboard the Spitzer 
Space Telescope, and compare them with 21 higher-luminosity quasar counter- 
parts. Silicate dust reprocessing dominates the mid-infrared spectra, and we 
generally measure the 10 and 18/im spectral features weakly in emission in these 
galaxies. The strengths of the two silicate features together are sensitive to the 
dust distribution. We present numerical radiative transfer calculations that dis- 
tinguish between clumpy and smooth geometries, which are applicable to any 
central heating source, including stars as well as AGNs. In the observations, we 
detect the obscuring "torus" of unified AGN schemes, modeling it as compact 
and clumpy. We also determine that star formation increases with AGN luminos- 
ity, although the proportion of the galaxies' bolometric luminosity attributable 
to stars decreases with AGN luminosity. 

Subject headings: dust — galaxies: active — galaxies: nuclei — infrared: galaxies 
— radiative transfer 



1. Introduction 

Accretion onto central supermassive (lO^-lO^^M©) black holes powers active galactic 
nuclei (AGNs). AGNs exhibit a great variety of observational characteristics, notably the 
presence or absence of spectrally broad emission lines, which determine their classification 



as type 1 or type 2, respectively. Unified AGN models (lAntonuccil 1 19931 ) account for these 
differences in terms of viewing geometry, with an optically and geometrically thick dusty 
torus that blocks the broad line region and central engine from some (type 2) lines of sight. 



The presence of the dusty torus can be detected not only in absorption but also in 
emission, with the bulk of the reprocessed AGN continuum emerging at infrared wave- 
lengths. The exact spectral energy distribution (SED) is a function of dust geometry, and 
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initial radiative transfer calculations model ed a smooth, uniform torus, which is consis- 



tent with basic requirements of uni fication (jPier &: Krolikl Il992l : iGranato &: Danesd Il994 



Efstathiou &: Rowan-Robinsorull995l ). However, rece nt high-resolution observations iiidicate 
that the obscuring torus has a clumpy geometry (jjaffe et al.l 12004 : iTristram et al.l 120071 : 
Beckert et al.ll2008l ). These data require a large range of dust temperatures to coexist at 
the same distance from the central heating source, whereas the temperature of a smooth 
torus declines with distance. New models that place the dust in a clumpy toroidal distri- 
bution account for these observations while remaining consisten t with unification schemes 
dNenkova et"alll2002l : iHonig et al.lbood : ISchartmann et aliboosh . 



Silicate dust specifically dominates the mid-infrared (MIR) spectra of galaxies, and it can 
produce both the continuum and prominent spectral features at 10 and 18/im. The stronger 
lO/xm feature originates from a SiO stre tching mode and the 18/im feature from a SiO 
bending mode (IKnacke fc Thomson! 1 19731 ). A view of a hot, optically thin surface produces 
the features in emission, and a view through a cold screen shows the features in absorption. 
A smooth torus therefore exhibits strong emission from type 1 orientations (which view 
the hot inner throat of the torus directly), and deep absorption in type 2 AGNs (viewed 
through the torus). However, observations of AGNs do not conform to these expectations. 
Silicate emission is not universally observed in t ype 1 galaxies, a n d only recently has it been 



detected, primarily in high-luminosity AGNs (iHao et al.l |2005| . 120071 : ISiebenmorgen et al. 



20051 : ISturm et al.ll2005l ). Current observations of s amples of the lo wer- luminosity Seyfert 1 



galaxies even show silicate absorptio n on average (IHao et al.l 12007), and Seyfert 2 galax ies 



typically show only weak absorption (IMason et al.l 12006 



Shi et al. 



20061 : iHao et al.ll2007f ). 



Here we use MIR spectroscopy of unobscured AGNs to diagnose their native dusty 
environments. Even more informative than the behavior of one silicate feature alone, the 
combination of both silicate features together reveals the geometry of the rep rocessing dust 



aroun d the AGNs, discriminating between smooth and clumpy distributions (ISirocky et al. 



20081 ). We also use the spectra to evaluate the energetic contribution of star formation, 
which is evident in low-ionization emission lines and polycyclic aromatic hydrocarbon (PAH) 
emission, and which accounts for a significant fraction of the long-wavelength continuum fiux. 



2. Sample and Data Reduction 



We select Seyfert 1 galaxies from the iRush et al.l (119931 ) 12/im survey, using 31 with 
archival low-re solution spectra from the Spitzer Space Telescope Infrared Spectrograph (IRS) 
in this study (IWerner et al.l l2004l : iHouck et al.ll2004l ). The hot dust continuum of AGNs 
dominates the MIR emission in this fiux-limited sample. We restrict the Seyfert Is to those 
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galaxies with optical classification between 1 and 1.5. As a comparison sample, we use 21 
nearby {z < 0.4) type 1 quasars having archival IRS spectra. We classify all sources with 
monochromatic luminosity {vLy) at 14/im L14 > 7 x 10^°Lq as quasars. We list the galaxies, 
their basic properties, and observational details in Table [TJ 

These IRS observations cover the MIR bandpass from 5.2 to 38/im. We used stan- 
dard Spitzer Science Center pipeline version S15.3.0 data and extracted the spectra with 
the Spitzer IRS Custom Extraction (SPICE) package. The observations were performed in 
either staring or mapping mode. We used two nodded positions of the same order to back- 
ground subtract the staring mode observations, in which the nucleus is always centered in 
the slit. Differencing exposures in first and second orders provided background subtraction 
of mapping mode data. For these mapping mode observations, in addition to the central 
spectrum obtained with the slit centered on the nucleus, we included contributions from the 
two adjacent off-center spectra. We compared the galaxies' FWHM to that of a calibration 
star to identify extended sources, finding eight: ESQ 12-G21, MCG -5-13-17, NGC 7469, IC 
4329A, NGC 1566, UGC 5101, NGC 3227, and MCG -6-30-15. We used default SPICE ex- 
traction for point sources and extracted only the central 4 pixels (which corresponds to 7.2" 
and 20.4" in the short- and long-wavelength orders, respectively) from extended spectra. We 
robustly averaged individual spectra from each order to remove bad data. We scaled short 
wavelength spectra to match the fiux of the long- wavelength order "Long Low 1," which has 
the widest slit and is less sensitive to pointing errors. 



3. Results 

3.1. Broadband Spectral Characteristics 

The Seyfert 1 sample and the comparison quasars reveal the typical MIR characteristics 
of unobscured AGNs alone, which are evident in their average spectra (Figures [T] and [5]). 
The individual spectra are normalized at 14/im, which represents the dust continuum, and 
weighted according to their signal-to-noise in each average spectrum. Both the Seyfert 
and quasar spectra show prominent lO/im silicate emission, which reveals the geometry of 
the dust distribution, high-ionization emission lines, which are predominantly the result of 
ionization by the AGN continuum, and polycyclic aromatic hydrocarbon (PAH) emission, 
which is associated with star formation. 

A few sample members are luminous or ultraluminous infrared galaxies (LIRGs or 
ULIRGs), having L/r > IO^^Lq or > 1 0^"^ Lc^, respectively. The se galaxies (plotted as 
red histograms) are preferentially mergers JSanders fc Mirabellll996h . and they often exhibit 
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extended MIR emission (jSoifer et al.ll2000l : lAlonso-Herrero et al.ll2006l ). which is not directly 
attributable to the AGN. We show the average Seyfert 1 and quasar spectra, both including 
(red) and excluding (black) these (U)LIRGs. The IR-luminous galaxies' spectra are dif- 
ferent from the others', showing silicate absorption and relatively strong long- wavelength 
emission. However, because only two of 31 Seyfert 1 galaxies and one of 21 quasars are also 
(U)LIRGs, the resulting average spectra are not significantly different in each case. We base 
the subsequent analysis on the (U)LIRG-free average spectra, which better isolate the AGN 
contribution. Table [2] contains the normalized average spectra of the Seyfert 1 galaxies and 
quasars. 

Figure [2] further shows both the average Seyfert spectrum (green) and the average 
quasar spectrum (black) together. Overall, these spectra are very similar in shape and 
emission features. The 5 to 14/^m flux density ratios, -F5/F14, are slightly but not significantly 
different, with F5/F14 = 0.25±0.12 in the Seyfert Is and -F5/F14 = 0.45±0.16 in the quasars. 
We can d escribe the MIR s pectra as a power law F„ oc u", where a < is characteristic 



of AGNs ( Elvis et al.lll994 



Klaas et al. 



2OOII : lAlonso-Herrero et al.ll2006l ). The power law 



indices range from -0.5 to -2.8 in the Seyfert 1 sample and -0.3 to -1.6 in the quasars, with 
a = —1.3 and —0.8 in the average Seyfert 1 and quasar spectra, respectively. These results 
agree with those of lAlonso-Herrero et al.l (120061 ) . who found —2.8 < a < —0.5 over the 
wavelength range of 3.6-8/im in AGN-dominated galaxies. 

We also measured the equivalent widths (EWs) of strong emission lines in the average 
spectra. We find the line EWs to be larger in the Seyfert galaxy spectra than in the quasar 
spectra, considering bo th strong high ex citation lines, such as [O IV] and [Ne V] , which are 
attributable to AGNs (ILutz et al.l 120031 ). and the low ionization lines [Ne II] and [S III], 
which originate in star formation. Specifically, the [O IV] EW is 0.09/im in the average 
Seyfert spectrum and 0.03/im in th e aver age quasar spectrum. These r es ults agree with 



previo us MIR work by iHonig et al.l (120081 ) and iKeremedjiev et al.l (120081 ) . iMelendez et al. 



( 120081 ) find L[oiv] oc -^^2-ioA;eV' which similarly indicates smaller EW in the higher luminosity 
AGNs, if we consider the 2-10 keV X-ray luminosity as a proxy for intrinsic AGN luminosity. 
Finally, while the quasar comparison sample was not selected r obust ly, we note the resulting 
average spectrum is extremely similar to that of iNetzer et al.l (120071 ) . showing a comparable 
spectral shape, broad silicate emission, and PAH emission. 



3.2. Type 1 AGN silicate emission 



Dust produces both the MIR continuum and the silicate features. We model the spec- 
tra assuming the same dust produces both, as opposed to invoking physically separate line- 
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and continuum-producing regions (c.f. ISchweitzer et al.ll2008l ). In the IRS spectra, we mea- 
sure the continuum over a short-wavelength region (typically 5-7/i.m), an intermediate point 
(around 14/im), and a long- wavelength region (typic ally 26.5-31.5 A tm), a nd we fit a spline 
to define the full continuum, utilizing the method of ISirocky et al.l (120081 ) . The silicate fea- 
tures are evident in the dust absorption cross section, and the local minimum around 14/im 
produces the pseudo-continuum of the observed spectra at this wavelength. The result- 
ing continuum fits agree well with the radiative trans fer model calculatio ns of the emission 



from synthetic dust that lacks the silicate features (ISirocky et al.ll2008l ). Figure [3] shows 



an example of the continuum fit to Mrk 766. The procedure slightly varies depending on 
the spectral characteristics, with "continuum-dominated spectra" containing AGN emis- 
sion lines and weak PAHs, "PAH-dominated spectra" exhibiting strong PAH emission, and 
"absorption-domiii a ted spectra" showin g strong ice and hydrocarbon absorption below 14/im 
(ISpoon et al.l l2005l : ISirocky et al.l 120081 ). Most of these spectra are continuum-dominated, 
which provide more reliable 10/im silicate measurements, while the silicate strength mea- 
surements of PAH-dominated and absorption-dominated spectra are more uncertain. 



We measure the silicate feature strength 



S, 



Sil 



In 



cont 



at the wavelength of the strength extremum around 10 and 18/im, w here Fphs is the observe d 
fiux density and Fcont is the fitted continuum fiux density, as in iLevenson et al.l (120071 ). 



Table El lists these silicate strengths (5*10 and ^ig) and peak wavelengths (Aio and Aig). The 
average values of Aio and Aig are 10.0 ± O.l/im and 18.1 ± 0.2yum in Seyfert 1 AGNs and 
10.1 ± 0.2jU m and 18.0 ± 0.2 / im iri quasars. These values are consistent with the dust cross 
sections of lOssenkopf et al.l (119921 ) and characteristic of the interstellar medium, although 
radiative transfer effec t s can result in small (< 0.3/im) wavelength shifts. In contrast, we 
note that ISturm et al.l (120051 ) attribute the large wavelength shifts they measure (Aio up to 
11.5/xm) to unusual grain size distributions. 

The silicate strengths reveal features in emission {Ssu > 0) in 49 of 52 AGNs in the 
combined sample (Figure Hj), although the emission is generally weak. The only type 1 
AGNs that show lO/xm silicate absorption are also LIRGs or ULIRGs. The corresponding 
average Seyfert 1 spectrum shows the 10p, m silicate feature in obvious emission (5*10 = 0.11), 
whereas the average Seyfert 1 spectrum of iHao et al.l (120071 ) shows weak silicate absorption. 
The homogeneous sample selection we employ and the exclusion of (U) LIRGs better isolates 
the AGN and its i mmediate environm ent in the MIR spectra. In contrast, the heterogeneous 
Seyfert 1 sample of iHao et all J2007h includes LIRGs and ULIRGs. Most of the red 2MASS 
sources, for example, are indeed LIRGs and ULIRGS based on their IR luminosities. The 
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dusty star- forming regions of these (U)LIRGs contribute significantly to their IRS spectra 
and alter the appearance of the silicate features. We conclude that the average Seyfert 1 
spectrum of Figure [H which shows lO/xm silicate emission, best characterizes the AGN and 
its immediate surroundings in the MIR. The higher luminosity quasars similarly show lO/xm 
emission, with 5*10 = 0. 18 on average. Sihcate emission in AGNs was first discover ed in 
high luminosity galaxies (IHao et al.ll2005l : ISiebenmorgen et al.ll2005l : ISturm et al.ll2005l ). and 
although the quasar sample shows stronger 10/im silicate strength than the Seyfert sample, 
we find no significant correlation of silicate strength with AGN luminosity. Finally, we find 
no trends with Aio or Ais, indicating that the chemical composition of the dust does not vary 
significantly with AGN luminosity or silicate strength. 



3.3. AGN-Star Formation Connection 



Indicators of star formation, including [Ne II] 12.8/im and several PAH bands (IGenzel et al 



19981 . and references therein), are present in the majority of the spectra. We measure the 
integrated luminosities of the [Ne II] and 6.2/im PAH emission to quantify the star forma- 
tion contribution, fitting a local continuum and a Gaussian for the line emission. The 5/im 
monochromatic continuum luminosity scales with the AGN luminosity, without the ambi- 
guity of a star formation contribution to the continuum luminosity that is present at longer 
wavelengths. These AGN and star formation luminosities are indeed positively correlated 
over both Seyfert 1 and quasar luminosities (Figure [5]). We certainly measure [Ne II] in 
all but two Seyfert galaxies, and the resulting robust [Ne Il]-5/xm correlation (dashed line) 
is consistent with the less complete quasar measurements. The [Ne II] non-detections are 
a consequence of poor signal-to-noise. Combining the spectra without detections, we suc- 
cessfully measure the line. We plot these sources (green) at this average [Ne II]-5/im ratio, 
considering the Seyfert 1 galaxies and quasars separately. The solid line shows the subse- 
quent correlation over all galaxies, including those in which [Ne II] is not detected directly, 
which agrees with the Seyfert 1 result using individual detections alone. 

The PAH measurements indicate similar trends of star formation increasing with AGN 
luminosity, but these results are less robust. We directly detect the PAH emission in fewer 
individual spectra (24 of 31 Seyfert Is and 8 of 21 quasars). While we recover the PAH 
emission in the average "PAH-less" Seyfert 1 spectrum, we do not certainly detect PAH 
emission in the average "PAH-less" quasar spectrum. Overall, this work agrees with that of 



Schweitzer et al.l (120061 ) who find a similar correlation between AGN and starburst luminos- 



ity, measuring 6/im continuum luminosity and PAHs in quasars. 



While luminosity due to star formation and accretion are correlated, the relative con- 
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tribution of star formation decreases with increasing AGN luminosity. The slopes of the 
linear fits are 0.50 ± 0.16 and 0. 66 ± 0.07 fo r the well-measured Seyferts and all sources, 
respectively. Using equation 1 of IHo &: Ketd (120071 ). we calculate the luminosity contribu- 
tion of the starburst component from the [Ne II] integrated luminosity, treating the infrared 
luminosity as an approximati on of bolom e tric lu minosity. Similarly, using the median SED 
and bolometric corrections of lElvis et al.l (119941 ). we obtain the bolometric scale factor for 
5//m continuum, L^oi = IO.IL5. For weaker AGN contributors (e.g., L5 = 10^'^Lq), the star 
formation luminosity is as much as 80% of the AGN contribution, whereas the star forma- 
tion luminosity of strong quasars (L5 = lO^^L©) is around 5% of their AGN luminosity. The 
contribution of star formation to the continuum flux increases with wavelength. Comparing 
the AGN-dominated 14/im flux density with that at 30/xm, at which dust-reprocessed stellar 
light becomes significant, we find the average F14/F30 = 0.44 ± 0.18 for the Seyfert Is and 
F14/F30 = 0.68 ± 0.27 for the quasars. Similar to the conclusions based on [Ne II] and PAH 
emission, these results also suggest that the relative luminosity of star formation is greater 
in the Seyfert 1 galaxies than in the quasars. 



4. Dust Geometry from Silicate Features 



The strengths of the 10 and 18/im silicate features together are sensitive to the distribu- 
tion of dust surrounding any heating source, including stars as well as AGNs. Smooth and 



clumpy distributions occupy distinct regions of the "feature-feature diagram" (ISirocky et al. 



20081 ). which shows Sis ^s. Siq (Figure [6]). We compare these radiative transfer calculations 
with the observations to discern the dusty environment of the AGNs we observed. While the 
two strength measurements in an individual galaxy cannot constrain all the free parameters 
of any of the models, the type 1 AGNs are located in an area of the diagram that only 
clumpy models occupy. 



Schweitzer et al.l (120081 ) alternatively model MIR spectra of AGNs with multiple inde- 



pendent emission components. Combinations of blackbodies represent the continuum, due 
to hot dust close to the nucleus. The more distant, cooler, optically-thin narrow line region, 
located at 100-200 dust sublimation radii produces the silicate features. A disadvantage of 
this approach is that it allows sili cate only in emi ssion. It never produces silicate absorption, 
which is typical of type 2 AGNs (IHao et al.ll2007l ). Indeed, assuming AGN unification, these 
models would instead pred ict silicate emissi o n from type 2 AGNs in general. Moreover, the 
hot inner edge of the disk ISchweitzer et al.l (120081 ) describe would also produce strong sili- 
cate emission that would be observed directly in unobs cured AGNs but i s abse nt from their 
model. The key difference between the description of ISchweitzer et al.l (120081 ) and ours is 
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the location of the sihcate-emitting region. Unfortunately, no current observations provide 
the spatial resolution to discriminate between them directly. Thus, we pursue here mod- 
els in which a common dust distribution simultaneously accounts for the MIR continuum 
and spectral features. This technique has the advantage that it may be consistently and 
directly applicable to obscured AGNs, in agreement with unified AGN schemes, although 
we acknowledge that additional separate emission regions may be present in some galactic 
centers. 

Independent of geometry, all optically thin configurations exhibit silicate emission and 
are located at the same point in the diagram, which corresponds to the feature strength in 
the optical cross section. Tracks of increasing optical depth move toward weaker emission. 
Only smooth dust distributions can exhibit large negative feature strengths, producing the 
temperature gradient that is essential for deep absorption. All smooth spherical distributions 
show similar slopes in the diagram, independent of th e density distributi on and total dust 
extent, which the dust's optical properties determine (ISirocky et al.ll2008l ). We plot several 
characteristic examples in Figure [61 Spherical distributions of clumps occupy a distinct 
region of the feature-feature diagram, never showing deep absorption, even for comparable 
total optical depth. The reason for this behavior is that both dark (absorbed) and bright 
(illuminated) cloud faces are observed in the clumpy distribution. The silicate eraission from 
the bright sides fills in the absorption trough, reducing its depth (INenkova et al.ll2002l ). 



We consider whether changes to the dust composition could provide smooth distributions 
that describe the observations on the feature-feature diagram. The dust's optical proper- 
ties determine both the locati on of the optically thin point and the slopes of the smooth 
model tracks on the d iagram (ISirocky et al.ll2008l ). C ompared to th e dust we employ here 
( Ossenkopf et al. 1992 ). the "astronomical silicate" of Draine ( 2003a| jbl). for example, shows 
a reduced 18/im feature relative to the lO/xm feature, so the silicate features of these smooth 
models (Figure 9 of Sirocky et al. 2008) are more similar to those of the observations. How- 
ever, the specific models that lie close to the data are optically thick, with ry ^ 10, which is 
inconsistent with the small optical depths measured in type 1 AGNs. This resulting optical 
thickness is a general problem for smooth models of any dust that exhibits the MIR silicate 
features. Only a contrived dust having extremely weak intrinsic silicate features could re- 
main optically thin in the observed region of the diagram, yet even such a forced solution 
could not then produce the observed range of strength ratios. We conclude that the data 
lie below the smooth model tracks for reasonable dust properties, a region that only clumpy 
models occupy. 



While the smooth spherical models offer a valuable contrast, they are inappropriate in 
these cases, failing to allow direct views of the central engine and the resulting spectrally 
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broad emission lines. Nevertheless, the smooth spherical geometries do usefully indicate 
some of the realm of the smooth torus of classical AGN unification schemes, which provide 
unobscured views along select lines of sight. Specifically, the obscured type 2 view through 
the torus is analogous to the spherical geometry, admitting no view of directly-illuminated 
hot dust. Strictly unobscured (type 1) lines of sight to the central engine view the silicate 
emission from the optically thin illuminated surface of the torus directly. The net result 
is silicate emission, located near the optically thin point of the feature-feature diagram. 



tribution from the cooler torus interior within the observing beam ( 


Granato & Danese 


1994; 


Efstathiou & Rowan-Robinson 


1995; 


Van Bemmel & DuUemond 


2003 


). In type 1 views, the 



sil icate strength i s sens itive to the shape of the inner edge of the dusty torus. The models 



of iPier &: KrolikI (119921 ). for example, do not typically exhibit strong silicate emission. How- 



ever, several independent lines of evidence argue against smooth dust distri butions generally . 



First, the MIR emission of AGNs is ob served to be effectively isotropic (ILutz et al.l 12004 



Buchanan et al.l l2006l : iHorst et al.l 120061 ) . whereas all smooth torus models predict type 1 
AGNs to be significantly brighter than type 2 AGNs in the MIR, for a fixed intrinsic lu- 
minosity. Second, all optically thick smooth des criptions produc e very deep absorption in 
obscured nuclei, which is not generally observed (iHao et al.l 120071 ). 



Thus, we pursue clumpy dust distributions, which iKrolik fc Begelmanl (119881 ) originally 
proposed. The clumpy geometry allows a large range of dust temperatures to coexist at 
the same distance, as opposed to the monotonic temperature decline with distance that 
is characteristic of smooth models. Interferometric observations of NGC 1068 with the 
VLTI, fo r example, resolv e the lO/im emission and indicate cool dust located close to the 
nucleus ( jJaffe et al.ll2004j ). Similarly, VLTI observation s of the Circinus galaxy and NGC 



3783 provide further evidence of a clumpy dust structure (ITristram et al.ll2007l ; iBeckert et al. 

mi- 



Initially, we minimize the number of model free parameters and consider a spherical 
distribution of clumps. Although this geometry does not generally provide any clear lines of 
sight to the central engine, the simplified spherical models are powerful, and they capture 
the essence of the MIR spectra. Fundamentally, the total population of clouds produces 
the observed MIR emission, and the SEDs are insensitive to viewing orientation effects 
in any case, even when the dust distribution is not spherically symmetric. (We explicitly 
demonstrate this result belo w, presenting calcul ations for toroidal distributions of clouds.) 
We foll ow the formalis r n of ISirockv et al.l (120081 ). which is based on the radiative transfer 
code of iNenkova et al.l (l2008al ). In the computations, the individual clouds are distributed 
according to Poisson statistics, with an average number of clouds along a radial ray, A^o- The 
clumps are radially distributed according to a power law, cx r~'^, from the dust sublimation 
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radius, Rd to an outer radius i?o, which we parameterize with Y = Ro/Rd- The bolometric 
luminosity of the central source sets Rd, with Rd - 0.4(Lto«/10^5erg ^-^f/'^ for the AGN 
heating spectrum and dust sublimation temperature of 1500 K. The optical depth of each 
cloud in the V band is ry, so the total average optical de pth through all c l ouds is NqTv- The 



dust includes both silicates and graphite, and we use the lOssenkopf et al.l (119921 ) cool silicate 
optical properties. The model results we present are applicable to any heating source, not 
restricted to AGNs, because the dust erases all signatures of the incident spectrum from the 
emergent MIR emission. 

We initially leave the model parameters unconstrained in order to show the effects 
different parameters have on the model curves in the feature-feature diagram. The plotted 
simulations therefore do not all represent best- (or even "good-" ) fitting models. In exploring 
these parameter variations, we will identify the parameter values that produce models that 
generally agree with the observations on the feature-feature diagram. 

Figure [6] shows the silicate strengths of the clumpy sphere models for a range of A^^o, 
with F = 30 and q = 1. The optical depth per cloud increases along each track of fixed 
A'o, from Ty = 10 (at the upper right) to 80. Overall, increasing Nq results in diminished 
silicate emission, with silicate absorption emerging for A^o > 3. Having more clouds increases 
the chance that bright faces are obscured, and views of dark, absorbed faces occupy more 
lines of sight. As the optical depth per cloud increases along the constant A"o track, silicate 
strength initially decreases. Around Ty = 60, which corresponds to tiq ^ 3, optical depth 
effects within individual clouds become important, and the lO/im silicate strength increases 
as the optical depth per cloud increases further. The 18/im strength generally continues to 
decrease, with rig < 2 per cloud when ry = 80. These clumpy models describe the type 1 
AGN data well, typically with a small number of clouds (A"o ~ 2). 

The MIR-emitting region is compact, and clouds located far from the AGN do not 
significantly affect the silicate feature strengths. We demonstrate this result first considering 
the radial density profile q = 2 models over a range of outer size, Y (Figure [7]). Having 
steep radial density profiles, these distributions are inherently compact and therefore are 
not sensitive to the outer extent. For example, 80% of clouds are located within 3.6Rd and 
4:.8Rd for y = 10 and Y = 100, respectively. However, the shallower radial distributions 
are sensitive to the total size because the number of nearby clouds is a function of both Y 
and Nq. In these cases, replicating the MIR behavior of the small- A'o compact distributions 
requires increasing Nq as the total size increases. 

Small numbers of clouds along radial rays in spherical clumpy models best match the 
silicate strengths of type 1 data, with A'o ~ 2. We find little variation in silicate strength as 
a function of other parameters when A^o is small, so the conclusion that the immediate sur- 
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roundings of these AGNs contain few clouds along radial rays is robust. However, for larger 
values of Nq, the silicate strengths depend sensitively on the combination of all parameters, 
including q and Y, as well. 

While the clumpy sphere captures the essence of the MIR emission from these AGNs, 
general support for unified AGN models, especially the requirement of unobscured hues 
of sight t o the central engine, f avors a toroidal distribution. We use the clumpy toroidal 
models of iNenkova et all J2002I . boOSal lbl). which allow variation of the viewing inclination 



angle, i, and torus scale height, a, in addition to the parameters of the spherical model. 
Figure [8] illustrates the clumpy torus. The quantity A'^o now represents the average number 
of clouds through an equatorial ray of the torus, and we consider Gaussian distributions, 
where the average number of clouds Nios{P) = Noexp{—i3'^/a'^) along angle P measured from 
the equator. (The inclination angle is measured from the symmetry axis of the torus, so 
(3 = 90° — i.) Unobscured views are more likely for small values of i, with the photon escape 
probability Pp„^ r. = e^vi—NioJP)) de scribing the likelihood of an unobscured view of the 



central engine (INenkova et al. 



2008tf ). 



The clumpy nature of the dust is fundamental, and as a result, the toroidal distributions 
occupy a region of the feature-feature diagram similar to that of the spherical arrangements. 
The toroidal calculations uphold the general conclusion that the MIR-emitting region is 
small. For example, these inherently-compact q = 2 results are insensitive to the total 
outer extent of the torus, similar to the spherical calculations above. The flat {q = 0) radial 
density profile, however, is a strong function of outer radius (Figure [H]). In this case, the most 
compact tori (having Y = 10) generally describe the data well, with A"o ~ 4. The clouds are 
spread over a large volume in the extended {Y = 100 or Y = 30) tori when q = 0. With this 
radial distribution few clouds are located close to the nucleus, even when the total number 
of clouds is large. For example, using model parameters Y = 100, q = 0, and A'o = 10, 
an average of only two clouds are located within 20/2^ along equatorial rays. As a result, 
this particular combination of parameters describes many of the observations. However, a 
significant fraction of observations lie above these models (having greater 5*18), even for large 
values of Nq. No change in other parameters can shift the q = 0, Y = 100 model tracks up 
to account for the stronger 5*18 measurements. Thus, this combination of parameter values 
is not generally characteristic of the observed sample, although it may describe particular 
galaxies. 

Computations of varying q for fixed Y again show that the nearby clouds determine the 
silicate feature strengths. In Figure [TDl we plot the strengths for models in which q ranges 
from to 2, with Y = 30,i = 30°, a = 45°, and A'o = 1, 2, 4, and 10. Models having steeper 
density distributions require fewer clouds along radial rays to match the silicate strength 
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of comparable models having shallower radial distributions. Furthermore, as q increases, a 
larger fraction of the clouds are located closer to the AGN. For example, with Y = 30, in 
a g = distribution, 80% of all clouds are located within 24Ra whereas the same fraction 
of clouds are confined to 4.4/?^ in the q = 2 density profile. The shallowest density profile 
{q = 0) does not follow the general trend of decreasing strength with increasing A'^o as 
rapidly as the steeper distributions do. Considering the observations, few clouds (< 6) along 
radial rays within a compact torus (size ~ ISi?^) accounts for the MIR silicate features of 
unobscured AGNs. 

Increasing the torus scale height reduces the silicate strength for a given Nq, as Figure [TT] 
shows. The total cloud distribution determines the behavior of the MIR emission, and an 
increase in a results in more clouds overall for a given A"o. Formally, Nq sets only the average 
number of clouds along equatorial rays, and a determines how rapidly the radial number 
declines with altitude. Again, with more clouds present, the bright cloud faces are more 
likely to be obscured, which reduces the net silicate emission and produces absorption in 
some cases. Here, the increased total number of clouds (Ntot) is a consequence of increasing 
a rather than increasing Nq alone. However, even for Nq = 10, the cr = 15° model shows 
emission, because N^ot is small in this narrow torus. Few high-altitude clouds are present to 
block views of the directly-heated cloud surfaces that exhibit silicate emission. 

The dominant direct view of a hot optically thin surface produces strong silicate emission 
in general, which is typical of smooth torus models of type 1 AGNs. The small-cx case 
approaches a two-dimensional dust distribution and explicitly demonstrates the failure of 
this simplification (in a clumpy or smooth arrangement), given the observed weak silicate 
features. Instead, the mixture of contributions from hot and cold cloud faces is the essence 
of the MIR emission. Thus, despite the defect of the spherical models in not providing 
unobscured lines of sight to the AGN in general, they better approximate the MIR results 
than a two-dimensional geometry does, and the development to the toroidal configuration 
ultimately corrects this problem. 

We consider several inclination angles for fixed Nq and a and find no significant differ- 
ences in the models' silicate strengths for small i typical of type 1 views {i < 40°). Thus, 
the silicate features do not usefully diagnose the viewing angle, and we adopt i = 30° in the 
comparisons below. According to unification schemes, the only difference between type 1 
and type 2 AGNs is the viewing angle, whereby direct views of the central engines of type 2 
AGNs are obscured, although the dusty AGN surroundings are inherently the same in both 
cases. A smooth torus strictly separates the different types at a particular viewing angle, 
distinguishing lines of sight through the dusty torus material from unobscured views near 
the symmetry axis. The clumpy formalism presents no strict dividing line, however. Instead, 
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type 1 or type 2 views may arise from any angle, but lines of sight near the equatorial plane 
are more likely to be obscured, and those near the symmetry axis are more likely to remain 
clear. The escape probability, Pesci describes the likelihood of a type 1 view, and it is a 
function of i, a, and Nq. 

We compare clumpy torus properties that describe the type 1 observations well to predict 
the silicate characteristics of their type 2 counterparts, viewed at higher inclination, which 
are obscured. We find that these type 2 AGNs show less silicate emission than type 1 views of 
the same dust distribution. Specifically, we calculate the sihcate strengths for combinations 
of model parameters that describe the type 1 measurements well: Nq = 2-4, Ty = 30-60, 
a = 30-60°, with Y = 10-30 for q = 1, and Y = 10-100 for q = 2. The corresponding type 
2 views {i = 70°) yield —0.4 < Siq < 0.15 and —0.15 < Sis < 0.18, a wider range of silicate 
strength than the i = 30° views exhibit. Overall, the silicate is generally weakly absorbed 
in the obscured AGNs, which agrees with previous observations ( iHao et al.l 120071 . Levenson 
et al. 2009, in preparation). Another way to test unified AGN schemes is to model the 
characteristic clumpy distribution of observed type 2 AGNs and compare this r esult with 



the ty pe 1 distributions. The peak of the lO/xm silicate strength distribution of the lHao et al. 



(120071 ) Seyfert 2s ranges over —0.4 < 5*10 < —0.1, which corresponds to A^o — 3-4 in the 



i = 70° clumpy torus models using the same parameter combinations as above. Changing 
inclination to i = 30° and holding all other parameters constant, the silicate strengths of the 
Aq = 3 or 4 models match those of the type 1 observations. These results are again consistent 
with standard AGN unification, whereby the central engine and its immediate environment 
are the same in all AGNs, and variations in viewing geometry alone produce observable 
differences. One further consequence of the clumpy geometry is that it c an produce silicate 



emission even when the AGN is obscured, which has been observed (ISturm et al.l 12006 



Teplitz et al.ll2006l : iHao et al.l 120071 ). unlike a smooth distribution, in which type 2 AGNs 
always exhibit silicate absorption. Even when the AGN and broad line region are hidden, 
direct view of some hot cloud faces can result in measurable silicate emission. 

Despite the model degeneracies, we identify several ranges of "standard" parameters 
that describe the data well and do not impose severe restrictions on other parameters. In 
particular, we favor Y = 30, q = 1, a = 45°, i = 30°, Ty = 30-60, and Nq < 6. Because the 
clouds close to the AGN determine the MIR emission, confining the total extent of the torus 
(Y = 10 or 30) generally produces models that describe the data well. Recent observa tions 



also show a cornpact MIR torus, with sizes of 3-5Rd and lO-R^, in iMason et al.l (120061 ) and 
Tristram et al.l (120071 ). respectively. The extended (Y = 100) torus is successful only when 
the clouds are concentrated (with radial density profile q = 2). Similarly, in the constant 
density distribution (g = 0), only Y = 10 yields acceptable results. 
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The MIR measurements alone do not strongly constrain the torus s cale height. Instead 



relative number s of type 1 and 2 Seyfert galaxies indicate a ~ 30 to 45° (ISchmitt et al.ll2001 



Hao et aljboosh. and consideration of the SED and lOfim silicate feature favors a 



30° 

(jNenkova et al.l l2008bl ) . However, the cr = 30° models underpredict Sis at lower values of 
5*10, even with large numbers of clouds {Nq > 10). Thus, we adopt a = 45° as the standard 
value and consider the range 30° < a < 60° to be applica ble to various individual galaxies. 



A disk-like geometry suggested by ISchweitzer et al.l (120081 ) with a clumpy distribution could 



be identified as a small-a torus, but low a values disagree with the observations. Because 
varying inclination angle does not significantly change the strength measurements provided 
that i < 40°, we discount directly face-on views and identify 20° < i < 40° to be typical of 
these type 1 AGNs. Furthermore, models of Ty = 30-60 produce silicate strengths that are 
similar to those of the observations. 

The physical characteristics of clumpy dust distributions are fundamentally different 
from those of smooth distributions. Both spherical and toroidal clumpy models cover the 
same regions of the feature-feature diagram, which remain inaccessible to all smooth descrip- 
tions. The behavior of the silicate features in the simplified clumpy spherical models yields 
results that are directly applicable to the more realistic clumpy torus models. To account 
for the MIR observations, the small-scale dust distribution is relevant, and both geometries 
show how the total torus extent and radial distribution together govern the effective com- 
pactness. In addition, small numbers of clouds within the small sphere or torus agree with 
the observations of type 1 AGNs. However, the Nq parameter is not constant across all 
models that describe the data: it increases from the sphere to the torus generally, and it 
increases with decreasing torus scale height. We therefore conclude that the total number 
of clouds available to reprocess the intrinsic AGN flux, Ntot (not Nq), and their distribution 
ultimately determine the behavior of the emergent MIR emission. 

Figure [12] demonstrates that for the same inputs of Y and q, fitting the data requires 
larger values of A'^o in the torus than in the sphere. In general, translating from any torus 
to sphere model requires NQ[torus) > NQ{sphere). The parameter Nq describes only the 
number of clouds along the equatorial ray of the torus, and the number of clouds diminishes 
with altitude. For a given Nq, the spherical model contains more clouds in total than the 
toroidal model does. Thus, in order to achieve the same N^ot in both geometries, Nq must be 
greater in the torus. The spherical models yield robust conclusions about the nature of the 
dust distribution around AGNs, confirming that few clouds within a small radius account 
for the observed MIR emission. However, because the spherical geometries generally fail to 
provide an unobscured view of the central engine, which these type 1 AGNs demand, we 
conclude that a clumpy, dusty torus characterizes the immediate surroundings of AGNs. 
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5. Conclusions 

Dust reprocesses the AGN continuum to emerge at MIR wavelengths, and we model 
both the resulting continuum and spectral features at 10 and 18/im due to a common dusty 
region. Isolating unobscured AGNs, we find these features in emission, both in Seyfert 1 
galaxies and in quasars. The emission is weak, however, with average emission strength 
SiQ = 0.11 and 0.18 in the Seyferts and quasars, respectively. In contrast, an optically thin 
medium, such as the directly-viewed hot interior of a smooth torus, would yield stronger 
emission (5*10 = 1.2). We conclude that the observable but weak emission is a consequence 
of clumpy AGN surroundings. We measure the peak wavelength Aio = 10.0 ± 0.1 and 
10.1 ± 0.2/im in the Seyfert Is and quasars, respectively. These values ar e consistent with 



radiative transfer computations using the optical properties of silicates that lOssenkopf et al. 



(I1992I ) model. In agreement with earlier work, we find that star formation increases with 
AGN luminosity, using [Ne II] and 6.2//m PAH to quantify the star formation contribution 
and the 5 and 14yum continua as proxies for the AGN luminosity. The fractional contribution 
of star formation to the total bolometric luminosity of these galactic centers decreases with 
AGN luminosity. Furthermore, considering the fiux ratio -F5/-F14, we can describe the type 1 
AGN spectra as a power law, with power law indices ranging from -0.5 to -2.8 in the Seyfert 
Is and -0.3 to -1.6 in the quasars. 

The strengths of the 10 and 18/im silicate features together are sensitive to the distribu- 
tion of dust surrounding any heating source, including stars as well as AGNs. We interpret 
the silicate strength measurements of these isolated AGNs as a consequence of a native AGN 
environment that is clumpy. The dust that determines the MIR behavior is confined to small 
scales. In radiative transfer calculations, either limiting the total radial extent of the dust 
or concentrating the cloud distribution effectively produces the compact distributions that 
describe the observations. 

A toroidal distribution of clouds is consistent with unified AGN schemes, offering un- 
obscured lines of sight to the central engine from some viewing angles. The total cloud 
distribution, not only the clouds located along the line of sight, determines the MIR emis- 
sion. As a consequen ce, the MIR outp ut is effectively isotropic, as observations of all types 



of AGNs show (e.g., iHorst et al.ll2006l ). We thus account for the model results in general: 
properties such as torus scale height and the number of clouds along radial rays that deter- 
mine the total cloud distribution govern the MIR behavior. Although spherical distributions 
do not generally provide the unobscured views these type 1 AGNs require, they do usefully 
capture the essence of the total cloud distribution while minimizing the number of free pa- 
rameters. While two silicate strength measurements cannot constrain all the free parameters 
of any of the models, we find that the type 1 AGNs are located in an area of the feature- 
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feature diagram that only clumpy models occupy. Few clouds are located along radial rays 
within the c ompact torus, which is consistent wi th column density variability observed in 
some AGNs (lElvis et al.l 12004 : iRisaliti et al.l 120071 ) . The MIR observations are not sensitive 
to the more distant cloud population, and we conclude the arrangement of dust immediately 
surrounding the AGN central engine is a clumpy torus that contains few clouds (< 6) along 
radial rays within a small radius (~ 15Rd)- 
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09 





.185 


4200 


11305728 


2005-APR-22 


480 


720 


PGl()48+342 


10 


51 


43. 


.9 


+33 


59 


27 


0. 


.167 


3800 


14192128 


2006-MAY-27 


480 


720 


PG1116+215 


11 


19 


08. 


.6 


+21 


19 


18 


0. 


.177 


4000 


4734464 


2004-MAY-14 


24 


56 


PG1121+422 


11 


24 


39. 


.2 


+42 


01 


45 


0. 


.226 


5300 


14193664 


2005-DEC-19 


480 


720 


3C273 


12 


29 


06. 


.7 


+02 


03 


09 


0. 


.168 


3600 


4978176 


2004-JAN-06 


168 


168 


Mrk 231 


12 


66 


14. 


.2 


+56 


62 


26 


0. 


.042 


880 


4978688 


2004-APR-14 


112 


120 


PG1307+085 


13 


09 


47. 


.0 


+08 


19 


49 


0. 


.155 


3500 


4735488 


2006-JAN-21 


24 


56 


PG13()9+355 


13 


12 


17. 


.8 


+35 


15 


21 





184 


4200 


4736000 


2004-MAY-14 


24 


56 


PG1322+659 


13 


23 


49. 


.5 


+65 


41 


48 





168 


3800 


14196224 


2006-,IUN-21 


480 


720 


IRAS F13349+2438 


13 


37 


18. 


.7 


+24 


23 


03 





108 


2400 


4373760 


2005-,TUN-07 


280 


280 


PG1352+183 


13 


54 


35. 


.6 


+ 18 


05 


17 


0. 


.152 


3400 


4736512 


2006-JAN-21 


24 


56 


PG1364+213 


13 


66 


32. 


.7 


+21 


03 


62 


0. 


.300 


7400 


14196992 


2006-FEB-02 


480 


720 


PG1402+261 


14 


06 


16. 


.2 


+25 


66 


36 


0. 


.164 


3700 


4676684 


2004-JUN-08 


66 


120 


PG1427+480 


14 


29 


43. 


.1 


+47 


47 


26 


0. 


.221 


5200 


14198528 


2006-JAN-18 


480 


720 


PG2130+099 


21 


32 


27. 


.8 


+ 10 


08 


20 


0. 


.063 


1400 


3761408 


2004-JUN-06 


168 


240 


PG2233+134 


22 


36 


07. 


.7 


+ 13 


43 


66 


0. 


.326 


8100 


4734208 


2003-DEC-17 


66 


120 



Note. — 



Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
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Table 2. Average Type 1 AGN Spectra 



Rest Wavelength 
(Mm) 
(1) 


Scyfcrt 1 Scaled Flux Density 

(Jy) 

(2) 


Quasar Scaled Flux Density 

(Jy) 

(3) 


5.000 


0.226 


0.441 


5.030 


0.237 


0.446 


5.060 


0.247 


0.449 


5.090 


0.257 


0.453 


5.120 


0.266 


0.450 


5.150 


0.271 


0.449 


5.180 


0.277 


0.453 


5.210 


0.287 


0.458 



Note. — Tablc[^is published in its entirety in the electronic edition of the Astrophys- 
ical Journal. A portion is shown here for guidance regarding its form and content. Col. 
(1): Rest wavelength. The wavelength scale is non-uniform and based on the orginal 
resolution, which decreases toward longer wavelengths. Col. (2) : Average Seyfert 1 
spectrum normalized at 14/j,m. Col. (3): Average quasar spectrum normalized at 14/im. 
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Table 3. Spectral Measurements 



Integrated Flux 

Aio Ais [Ne II] 6.2^m PAH F5 F14 F30 

Galaxy Sio (m™) Sis (pm) (erg s~'^ cm~^) (erg s"-"^ cm~^) (Jy) (Jy) (Jy) 



Seyfert Is 



Q' 



ESO 12-G21 


0, 


.15 


10. 


,1 


0, 


.08 


18. 


,2 


2, 


lE- 


-13 


5. 


9E- 


■ 13 


0. 


,05 


0, 


.15 


0-34 


Mrk 335 


0, 


,15 


10. 


,0 


0, 


10 


18 


,3 














0. 


,11 


0. 


22 


0.35 


iVirK ±Uo^ 






























0. 


,02 








NGC 931 





.05 


10 


,1 


0, 


08 


18 


,2 


8. 


6E- 


-14 


9. 


6E- 


-14 


0. 


,16 


0, 


54 


1.10 


IRAS F03450+0055 


0. 


.14 


9 


.9 


0. 


06 


18 


,0 














0, 


,09 


0. 


31 


0.49 


NGC 1566 


0. 


.02 


10. 


.0 


0. 


16 


18. 


2 


2. 


2E- 


-13 


Y, 


6E- 


-13 


0, 


.06 


0. 


15 


0.43 


3C120 


0. 


.26 


10. 


,2 


0. 


15 


17. 


.8 


6. 


3E- 


-14 


1, 


OE- 


-13 


0, 


.09 


0. 


31 


0.63 


MCG -5-13-17 





.09 


9 


.8 


0. 


10 


18 


,2 


1, 


6E- 


-13 


2. 


3E- 


-13 


0, 


.03 


0. 


21 


0-64 


Mrk 6 





.24 


10 


2 


0, 


18 


18 


I 


2, 


4E- 


■13 


2. 


OE- 


-13 


0, 


.11 


0, 


32 


0.70 


Mrk 9 


0. 


.04 


9 


,7 


0. 


14 


17. 


.8 


7, 


8E- 


■ 14 








0, 


.09 


0. 


24 


0.52 


Mrk 79 


0. 


.07 


10. 


,1 


0. 


10 


18. 


.0 


7, 


OE- 


-14 








0, 


.15 


0. 


48 


1.07 


Mrk 704 


0. 


.07 


10. 


,2 


0. 


05 


18. 


,2 


2. 


9E- 


-15 








0, 


.16 


0. 


44 


0.47 


UGC 5101 


_1, 


.52 


9 


.9 


— 0. 


.19 


18. 


.3 


6. 


.4E- 


-13 


1, 


4E- 


-13 


0, 


.09 


0. 


.27 


2.24 


Mrk 1239 


0. 


.13 


10. 


,1 


0. 


11 


18. 


.3 


1, 


5E- 


-14 


2. 


5E- 


-13 


0, 


.36 


0. 


82 


1.37 


NGC 3227 





.01 


10 


.0 


0, 


06 


18. 


.2 


1, 


lE- 


■ 12 


1. 


6E- 


■ 12 


0, 


.13 


0, 


63 


1.91 


NGC 3511 





.07 


1(J 


2 


0, 


07 


18. 


,2 


3, 


2E- 


-13 


3. 


8E- 


-13 


0.004 


0, 


()6 


(J. 22 


NGC 3516 


0. 


.03 


10. 


.0 


0. 


10 


17. 


,7 


9. 


9E- 


■15 


9. 


3E- 


■ 14 


0, 


.12 


0. 


39 


0.92 


NGC 4051 


0. 


.07 


10. 


.0 


0. 


.06 


18. 


.0 


3. 


.lE- 


-13 


7, 


2E- 


-13 


0. 


.14 


0. 


.64 


1.43 


NGC 4151 


0. 


.14 


10. 


.0 


0. 


.15 


18. 


4 


1, 


3E- 


-12 


6. 


2E- 


-13 


0, 


.66 


2. 


.45 


4.06 


Mrk 766 


0. 


.05 


10. 


,1 


0. 


08 


18. 


.0 


3. 


.3E- 


-13 


3. 


6E- 


-13 


0, 


.11 


0. 


.57 


1.75 


NGC 4593 


0. 


.11 


10. 


.0 


0. 


08 


18. 


,1 


6. 


8E- 


-14 


1, 


6E- 


-13 


0, 


.16 


0. 


46 


0.97 


MCG -2-33-34 





.04 


9 


.8 


0, 


15 


17 


.8 


1, 


2E- 


■13 


2. 


lE- 


-13 


0, 


.02 


0. 


12 


0.39 


MCG -6-30-15 





.02 


10 


,2 





10 


18 


3 


1, 


7E- 


-14 


1, 


OE- 


-13 


0, 


.13 


0, 


45 


0.79 


IC 4329A 


0. 


.04 


9 


.9 


0. 


04 


18. 


.1 


2. 


.4E- 


■13 








0, 


.25 


1. 


36 


2.05 


NGC 5548 


0, 


.14 


10, 


.1 


0. 


.08 


18. 


.1 


1. 


.2E- 


-13 


7. 


.9E- 


-14 


0, 


.06 


0. 


.29 


0.69 


Mrk 817 


0, 


.08 


10, 


.1 


0. 


.07 


18, 


.0 


2. 


.4E- 


-14 


1. 


.6E- 


-13 


0, 


.10 


0. 


.40 


1.36 


NGC 6860 


0, 


.06 


10, 


.1 


0. 


.10 


18, 


.0 


6. 


.9E- 


-14 








0, 


.11 


0. 


.25 


0.36 
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0, 


.12 


10, 


.1 


0. 


19 
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.0 


1. 


.6E- 
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3. 
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0, 
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0. 


.35 


0.66 
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0, 
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.1 


0. 


16 
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.9 


2. 


8E- 
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7. 
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0, 
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0, 
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0.46 
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9, 


.8 


0, 


11 
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8 


2, 
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4. 
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1. 
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7.96 
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0. 


.01 


17, 


.8 


2. 


.3E- 


• 14 


8. 


6E- 


■ 14 


0, 


.16 


0. 


,30 


0.84 


PG0947-H396 


0, 


.09 


10, 


.3 


0. 
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Fig. 1. — Average Seyfert 1 spectra (smooth heavy hnes) show 10 and 18/im sihcate emission. 
The bold black line best isolates the AGNs alone and is the average of 29 individual sources 
(grey histograms, normalized at 14/im) excluding LlRGs UGC 5101 and Mrk 1034 (red 
histograms). The bold red line is the average spectrum of all sources. In addition to the 
silicates, emission features that are characteristic of both AGNs and star formation are 
labeled. 
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Rest Wavelength [|im] 



Fig. 2. — Average quasar spectra (smooth heavy hnes) show 10 and 18/im sihcate emission. 
The spectrum plotted in black best isolates the MIR emission of these luminous AGNs and 
excludes the ULIRG Mrk 231 (red histogram), which contributes to the average plotted in 
red. As a comparison, the average Seyfert 1 spectrum is overplotted in green. The equivalent 
widths of both AGN-originating emission lines, such as [O IV], and star formation indicators, 
such as PAHs and [Ne II], are larger in the average Seyfert 1 spectrum. 
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Fig. 3. — Spectrum of the typical Seyfert 1 galaxy Mrk 766 illustrates the continuum fitting 
technique. The resulting continuum (dotted red) is a spline fit to measurements over a 
short-wavelength region, an intermediate point around 14/im, and a long-wavelength region 
(solid red). Vertical lines mark the measured peak wavelengths of the 10 and 18/im silicate 
features, which both appear in emission here, although the 10/xm emission is very weak. 
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Fig. 4. — Feature strengths at 10 and 18/im (S'lo and S\^. Positive strengths show sihcate 
emission in 49 of the 52 type 1 AGNs. The only type 1 AGNs that do not show 10/xm sihcate 
emission are also LIRGs or ULIRGs. Vertical and horizontal lines at zero strength separate 
the regions of emission [S > 0) and absorption [S < 0). Triangles identify Seyfert 1 AGNs, 
and circles mark quasars. 
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Fig. 5. — AGN and star formation luminosities, which the 5/im monochromatic continuum 
luminosity and the [Ne II] 12.8/im integrated luminosity indicate, are positively correlated, 
although star formation becomes proportionally less important with increasing AGN lumi- 
nosity. We certainly measure [Ne II] in all but two Seyfert Is and six quasars. Considering 
the Seyferts and quasars separately and averaging spectra without detections, we estimate 
[Ne II] in these cases (green). The solid line shows the correlation over all galaxies, which 
agrees with the Seyfert 1 result using certain measurements alone (dashed). Other symbols 
as in Figure m 
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Fig. 6. — Together, 10 and 18/im silicate feature strengths are sensitive to the dust geometry, 
and the AGN observations reveal clumpy surroundings. Independent of geometry, all opti- 
cally thin configurations are located at the same point in the diagram (square), and tracks 
of increasing optical depth move toward weaker emission and eventually show absorption. 
All tracks of smooth spherical dust distributions (black and red), have similar slopes and 
lie separate from the data, for a range of spatial extent, Y, and radial density distribution 
(oc r"*?). Models of clumpy environments (green) occupy regions of this "feature- feature di- 
agram" that are inaccessible to the smooth models. Each clumpy sphere track is a function 
of the average number of clouds along radial rays, A^O) with q = I and Y = 30. Clumpy 
models of Aq ~ 2 generally agree with the data. 
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Fig. 7. — Clumpy spherical distributions as a function of radial extent, Y, for q = 2. For 
Y = 10, 30, and 100, we plot curves of Nq = 1, 2, 4, and 10. Because the clouds closest to 
the AGN determine the MIR behavior, this inherently compact steep {q = 2) distribution is 
insensitive to the outer radius. 
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Fig. 8. — Cartoon of the clumpy torus. Emission and obscuration of the central engine are 
functions of viewing angle, i. Clouds are distributed from the dust sublimation radius, Rd, 
to the outer radius, Ro, according to a radial power law. The average number of clouds along 
an equatorial ray is Nq. The scale height of the distribution is a, with the average number of 
clouds Nios{P) = NQexp{—i3^/a^) along a radial ray at angle /3 measured from the equator. 
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Fig. 9. — Clumpy torus models as a function of Y , for q = 0. For Y = 10, 30, and 100, we 
plot curves of A^^o = 1, 2, 4, and 10, fixing a = 45° and i = 30°. In this distribution that 
is constant with radius, the size of the torus and Nq are related, with increasing values of 
both Y and A^o together producing results similar to models having smaller size and cloud 
number. The clouds close to the AGN govern the MIR emission, so in a large torus, A^o must 
increase to provide enough clouds at small radius to match the features of a smaller torus 
having fewer clouds that are all confined to the small scale. 
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Fig. 10. — Clumpy torus models as a function of radial density profile, q. We plot curves 
for A^o = 1; 2, 4, and 10, fixing F = 30 and i = 30°. Increasing q places more clouds close 
to the AGN for a given value of A^o- Thus, the models having steeper density distributions 
require fewer clouds to match the silicate emission of the bulk of the data. The shallowest 
density profile (g = 0) does not follow the general trend of decreasing strength as rapidly as 
the steeper distributions with increasing A"o. A compact torus (effective size ~ ISi?^) having 
few clouds along radial rays (< 6 along the equator) accounts for the MIR silicate features 
of unobscured AGNs. 
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Fig. 11. — Clumpy torus models as a function of a. We plot curves of A^^o = 1, 4, and 10, 
fixing Y = 30, q = 1, and i = 30°. Increasing the torus scale height reduces the silicate 
strength for a given A^^o, obeying the general trend of decreasing strength with increasing 
total cloud number. However, even for A^^o = 10, the cr = 15° model shows emission, because 
this thin torus contains fewer clouds to block direct views of bright cloud faces. We measure 
the same effect for other values of Y and q. 
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Fig. 12. — Clumpy torus and clumpy sphere models. We plot curves of A^o = 1, 2, 4, and 
10, fixing y = 30 and g = 1 in both cases, and i=30° and a = 45° in the torus models. The 
sphere and torus results are extremely similar for different values of Nq, with larger Nq in 
the torus. Because all clouds reprocess the AGN light and contribute to the MIR emission, 
the larger value of the parameter Nq in the toroidal distribution is required to have the same 
total number of clouds as the spherical distribution with smaller Nq. 



